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In recent years we have witnessed a growing theoretical interest to understand the 
critical properties of non-periodic (quasi-periodic, aperiodic or random) systems [1-4]. Pe- 
riodic systems, which are built of regularly spaced blocks or unit cells, are assumed to 
share the critical properties of the homogeneous model. Indeed close to the critical point 
the correlation length is much larger than the linear size of the blocks, thus averaging over 
correlated domains the periodic system becomes homogeneous. This argument naturally 
fails, if the size of the unit cell is diverging, i.e. for non-periodic models. In these systems 
- in the spirit of the Harris criterion[5] - one should investigate the fluctuations of the cou- 
plings in a domain of correlated spins. The perturbation, caused by the non-periodicity 
is then irrelevant (relevant), if the local energy fluctuations are smaller (greater) than the 
corresponding thermal energy. 

In study the critical behaviour of specific systems it is very important to know the 
position of the critical point exactly, since with this information the accuracy of the nu- 
merical methods is largely increased. In two-dimensions the position of the critical point 
of some homogeneous models can be deduced from duality relations[6]. The same method 
works for a random bond Ising model on the square lattice, where the two types of bonds 
are distributed randomly with equal probability [7], [8]. For periodic Ising models on the 
square lattice with unit cells of size m x n an exact relation for the critical temperature 
can be derived with the use of the Pfafnan method[9] or the transfer matrix technique[10]. 
This relation is given by a matrix equation which involves m-product of non-commuting 
2 n ~ 1 x 2 n_1 matrices. Therefore in the non-periodic limit m — >■ oo the critical temperature 
is explicitly known only for n = 1, i.e. for simple layered systemsfll]. 

The homogeneous Ising model has also been solved exactly on triangular and hon- 
eycomb lattices for a long time[12],[13]. The partition function of the models on these 
hexagonal lattices are related through a star-triangle transformation [6]. A repeated use of 
this mapping is the basis of an exact renormalisation group transformation by Hilhorst, 
Schick and van Leeuwen [14], which method has also been used to study the surface critical 
behaviour of layered models[15]. With this method one can also obtain numerical estimates 
on the critical point of layered Ising models on hexagonal lattices, however up to now the 
exact position of the phase transition temperature is not known analytically. 

In the present paper we consider this problem and derive an explicit expression for the 
critical temperature for arbitrary distribution of the couplings both on the triangular and 
honeycomb lattices. In the calculation one may directly generalize Houtappel's solution 
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for the homogeneous model[12]. We will, however, proceed on a simpler way and utilize 
the correspondences between the square and hexagonal lattices. Let us consider a square 
lattice, where every second horizontal couplings have a special strength, denoted by K 
(Fig la). Taking the limit K — > oo the lattice becomes equivalent to the triangular lattice, 
whereas for K = one is left with a honeycomb lattice. Thus the exact solution of layered 
triangular and honeycomb lattice Ising models can be obtained from the corresponding 
results of periodic models with a m x 2 unit cell on the square lattice. 

The problem of periodic Ising models on the square lattice has been solved analytically 
first by Hamm[9] with the Pfaffian method and later by Hoever[10] with the transfer matrix 
technique. Here we use results and notations of the second method. 

Let us consider M x N Ising spins on a square lattice. As indicated on Fig. la the 
vertical couplings Kj (j = 1,2,... M, k = 1, 2, . . . N) are given by: 



The model is periodic in the horizontal direction: Kj +m = Kj, Kj +m = Kj. In the 
following we consider only purely ferromagnetic models, i.e. Kj > 0, Kj > 0. As men- 
tioned before K — > oo and K = correspond to the triangular and honeycomb lattices, 
respectively. Introducing the notations 




(la) 



whereas the horizontal couplings are the same within one column: 




(16) 



Xj = sinh(2Kj) = (sinh(2Lj)) 



-l 



(2) 



one can write the partition sum of the model as [10], [11]: 



m 



n iV M 
(2xj) 2m trace T 



(3) 



where T is the transfer matrix: 



m 



T=l[T(j), T(j)=T 1 (j)T 2 (j) 



(4) 



Here 




(5) 
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whereas T 2 (j) is different for odd and even indices: 

/n/2 N/2 

T 2 (2Z - 1) = exp ( K2i-i<T Z 2k-i°2k + ^Vffc+i 
,fc=i k=i 

(6) 

N/2 N/2 

T 2 (2Z) = exp ( K^lk-i^k + J2 K ^ a 2k+i 
\ k=i k=i 

and cr£, a\ are Paufi matrices at site k. 

After performing the Jordan- Wigner transformation the transfer matrix is expressed 
in terms of fermion operators, which in Fourier representation reads as: 

m 

T= II T ^ T (Q) = I[T j ( q ), T(j)= J] Tj(q) (7) 

0<q<n/2 j = l 0<g<7r/2 

where the allowed values of q are equidistantly spaced by 2n/M. Denoting the leading 
eigenvalue of the transfer operator T(q) by exp(A(g)), the free energy per spin / is given 
by: 

/2 

(8) 



_ /3/= _x:,„ ( 2 % ) + — j o dq A (q) 



r- 

At the phase transition temperature the free energy is singular, which is a consequence 
of the singularity in X(q) at the border of the Brilliouin-zone, i.e. at q c = 0. Our aim in 
the following is to determine the position of this singularity, i.e. to calculate the critical 
temperature of the system. 

First, following Ref[10] we introduce the parity operator 

P = ofo% (9) 

which commute with the transfer operator T(g), thus T p (q) acts only on a subspace with 
parity eigenvalue P = ±1. For sufficiently high temperatures the largest eigenvalue A(g c ) 
is in the sector with P = 1, whereas for T = it is usually in the sector with P = — l.f 
Since the jump in the largest eigenvalue from one sector to the other can not be done 
analytically there is a phase transition point in the system at a temperature T c where 

\ + (q c ,T c ) = \-(q c ,T c ) (10) 



| It could happen in frustrated systems with ferro- and antiferromagnetic couplings, 
that at any finite temperatures X + (q c ) > ^~(q c )- Then the free energy is analytical and 
there is no phase transition in the system. 
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At the border of the Brillouin-zone the transfer operator Tj (g c ) can be easily evaluated. 
It is a 2 x 2 matrix given by [10]: 

where C 3 = cosh[Lj(l + P)], Sj = sinh[Lj(l + P)] and Ej = exp [(K - K J P)(-Py] . Thus 
the transfer operator is different for odd and even rows. 

In the following we consider the triangular lattice, which corresponds to K — > oo. In 
this limit the transfer operator has a simple structure. Keeping terms in leading exponen- 
tial order the product of two consequtive transfer operators is given by: 

nV(9c)= 

j'=i 

and 

n^)= 

j'=j 

for the P = 1 and P = — 1 sectors, respectively. Similarly, in this limit any product of the 
Tf (<? c )-s is given in a simple form and the eigenvalues X ± (q c ) can be expressed explicitly. 
Then from eq(10) one obtains the criticality condition as: 

m 

Y[ sinh(2Kj) exp(2Kj) = 1 (13) 

This remarkably simple formula is similar to the critical point condition obtained on the 
square lattice for a simple layered model, i.e. for n = 1 [11] . It is easy to check that eq(13) 
contains, as a special case, the criticality condition for the homogeneous system[12]. 

Next we turn to study the critical point condition for the honeycomb lattice. Let us 
first consider that system which is obtained from the square lattice on Fig la by removing 
the wavy bonds. The transfer operator Tj (q c ) of the system, given in eq(ll) with K = 0, 
is still a 2 x 2 matrix, i.e. not essentially simpler than for the original problem on the 
square lattice. Therefor there is no explicit relation, similar to eq(13), for the critical point 
of the layered honeycomb lattice Ising model which corresponds to Fig la with vanishing 
K bonds. However, if the layered structure is that shown on Fig lc, i.e. the vanishing 
bonds on the corresponding square lattice are horizontal, then one can derive a similar 
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] J exp(K) exp(— ifj/)sinh(2Lj/) 



coth(2L j+ i) 1 



+ 0(1) (12a) 



i+i 

jf[ exp(K) exp(Kj') 



1 




+ 0(1) 



(126) 



relation to eq(13). Here, instead of repeating steps of the previous calculation we make 
use the star-triangle mapping[6] and transform the result in eq(13). This transformation 
maps the honeycomb lattice with bonds pj, pj onto a triangular one with lattice constants 
Kj and Kj, such that 

exv(-4K )~ C ° sh2(p) cxpf ij^ ) - cosh (Pj + 2 Pj) (u) 

exp{ 4 ^J- cosh(pj+2 p. )cosh(pj _ 2 p j) > ex Pl 4A ^ " CO sh( Pj - 2pj) [M) 

Substituting these relations into eq(13) one obtains the criticality condition for the layered 
honeycomb lattice: 

m 

[J sinh(2pj)tanh(pj) = 1 (15) 

i=i 

As a special case this equation contains the criticality condition for the homogeneous Ising 
model on the honeycomb lattice[12]. 

We note that relations in eqs(13) and (15) have been numerically verified by the re- 
cursion method of Hilhorst and van Leeuwen[15]. Results on the surface critical behaviour 
of aperiodic Ising models, obtained by that method will be presented in a separate publi- 
cation^]. 
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Figure Captions: 



Fig 1 Connection between the square (a) triangular (b) and honeycomb (c) lattices. The 
square lattice with infinitely strong K bonds (denoted by wavy lines) is equivalent to 
the triangular lattice, which is related to the honeycomb lattice through a star-triangle 
transformation^] . 
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